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Transverse momentum correlation in azimuthal angle of produced hadrons due to mini-jets are 
studied first within the H1JING Monte Carlo model in high-energy heavy-ion collisions. Jet quench- 
ing in the early stage of thermalization is shown to lead to significant diffusion (broadening) of 
the correlation. Evolution of the transverse momentum density fluctuation that gives rise to such 
correlation in azimuthal angle in the later stage of heavy-ion collisions is further investigated within 
a linearized diffusion-like equation and is shown to be determined by the shear viscosity of the evolv- 
ing dense matter. Such a diffusion equation for the transverse momentum fluctuation is solved with 
initial values given by HIJING and together with the hydrodynamic equation for the bulk medium. 
The final transverse momentum correlation in azimuthal angle is calculated along the freeze-out 
hyper-surface and is found further diffused for larger values of shear viscosity to entropy density 
ratio rj/s ~ 0.2 — 0.4. Therefore the final transverse momentum correlation in azimuthal angle can 
be used to study the thermalization of mini-jets in the early stage of heavy-ion collisions and the 
viscous effect in the hydrodynamic evolution of the strongly coupled quark gluon plasma. 



1. Introduction. In high-energy heavy- ion collisions 
at the Relativistic Heavy-ion Collider (RHICV enormous 
amount of transverse energy is produced 0, i, Q via 
hard and semi-hard collisions of beam partons from each 
colliding nuclei [J, HI, d, 0, Final state interaction 
among these produced partons during the early stage of 
heavy-ion collisions will lead to a thermalized (or par- 
tially thermalized) quark-gluon plasma (QGP) which will 
undergo further hydrodynamic evolution and expansion 
before hadronization and freeze-out. The strong ellip- 
tic flow as measured in non-central Au + Au collisions 
at the RHIC energy Q indicates indeed the early ther- 
malization of the dense partonic system [Io| before hy- 
drodynamic expansion takes place and the momentum 
anisotropy of the medium partons is developed. The 
phenomenal jet quenching as observed in the suppres- 
sion of large transverse momentum single hadron spectra 
pd| and back-to-back hadron correlation [l2j also points 
to strong interaction between energetic partons and the 
bulk medium in the early stage of collisions. Extrapola- 
tion of the jet quenching picture to low and intermediate 
energy partons leads naturally to a picture of strong par- 
ton interaction in the bulk medium, a prerequisite for 
parton thermalization during the early stage of the for- 
mation of QGP in heavy-ion collisions. 

Jets in high-energy collisions are characterized exper- 
imentally by a collimated distribution of hadrons in the 
direction of the underlying energetic parton and are mea- 
sured as such in the calorimetric study of final transverse 
energy distribution. Such a calorimetric study of jets be- 
comes problematic and eventually impossible when the 
transverse energy of jets becomes increasingly small and 
comparable to the background of the underlying soft in- 
teraction. However, the characteristic particle correla- 



tion within a jet and back-to-back correlation from di-jets 
should still be present in the azimuthal angle correlation 
of produced hadrons if mini-jets dominate the hadron 
production mechanism (l3| . Such correlation due to jet 
and mini-jet production in heavy-ion collisions has been 
used to study jet quenching 1 121 and possible medium 
response to jet propagation [l4[ at RHIC. Extension of 
such study of correlation to the bulk matter should shed 
light on mini-jet thermalization throughout the evolu- 
tion of the dense matter. Since correlations in a ther- 
mally equilibrated medium also carry information about 
the dynamics in an interacting QGP [ill [H, [T3|, study 
of the transverse momentum correlation in relativistic 
heavy ion collisions in principle can also provide con- 
straints on the shear viscosity of the interacting dense 
medium [lU, 0, |2(j. Furthermore, experimental study of 
back-to-back hadron correlation in azimuthal angle can 
also provide information on the initial parton produc- 
tion mechanism and whether mono-jet production in the 
gluon saturation model [21| dominates. 

In this paper, we propose and study the transverse mo- 
mentum correlation in azimuthal angle as a direct mea- 
sure of both the early parton thermalization and later 
viscous hydrodynamics evolution in high-energy heavy- 
ion collisions. We first calculate the transverse momen- 
tum correlation in azimuthal angle within the HIJING 
[221 ] Monte Carlo model and study the effect of early par- 
ton interaction through jet quenching mechanism. We 
will further study the evolution of transverse momentum 
fluctuation and correlation in azimuthal angle within the 
framework of viscous hydrodynamics and the diffusion of 
the transverse momentum correlation due to shear vis- 
cosity of the bulk matter. 

2. Transverse momentum correlation. Instead of 
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hadron multiplicity correlation, we will focus on the 
transverse momentum correlation for the purpose of our 
study, which should be less sensitive to the hadronization 
process and can be directly related to the momentum 
density fluctuation in an diffusion equation derived from 
viscous hydrodynamics [3 13 as we will show later. We 
integrate over a finite size of psudo-rapidity \rj\ < 1.0 and 
discretize the azimuthal angle (f> into a number of bins 
whose size 8(f) <C 1 should be much smaller than the typ- 
ical size of a jet. The transverse momentum correlation 
between two bins centered at azimuthal angle 4>\ and (f>2 
is defined as 



(SEtiSEt2) 

(Eti) {Et2) 



(1) 



where SErj = Exj — (Erj) (j — 1,2) are trans- 
verse momentum fluctuations in the j-th bin, Etj — 
*^2iPTi8{8(f)/2 — \4>i — <t>j\) ar e the total transverse mo- 
menta in each bin, and the average is taken over events. 

We will first study the above transverse momentum 
correlation using the HIJING Monte Carlo model [22T ]. 
which is based on a two-component description of soft 
and semi-hard interactions. Mini-jets contribute to a sig- 
nificant fraction of hadron production at the RHIC en- 
ergy. Therefore, transverse momentum correlation Ct, 
as defined in Eq. (JTJ) , should reflect that of hadrons from 
mini-jets as shown in Fig. [TJ as the dashed line for cen- 
tral Au + Au collisions at the RHIC energy without jet 
quenching. The transverse momentum correlation is sim- 
ilar to two-hadron correlation with finite transverse 
momentum cut-off. It contains both near-side correlation 
due to hadron production from a single mini-jet and back- 
side correlation due to hadron production from two back- 
to-back mini-jets. Final-state interaction in HIJING is 
modeled through a schematic parton energy loss mecha- 
nism, which will lead to a broadening or diffusion of the 
transverse momentum correlation as shown by the solid 
line in Fig. [TJ Since the parton number density decreases 
very fast due to longitudinal expansion, mini-jet inter- 
action is the strongest in the early stage and therefore 
jet quenching represents an early thermalization process 
for mini-jets in heavy-ion collisions. This will lead to 
significant modification of the characteristic correlation 
due to mini-jets. Therefore, the transverse momentum 
correlation in heavy-ion collisions can be used as a direct 
measure of mini-jet thermalization. Since the HIJING 
results without jet quenching is equivalent to that of a 
superposition of independent p + p collisions, one there- 
fore should compare heavy- ion collisions to p+p to study 
the degree of mini-jet thermalization. As illustrated by 
the HIJING results in Fig. [1] early parton interaction 
leads to a significant diffusion of the near-side correla- 
tion while the shape of the away-side correlation remains 
almost intact since it is mainly controlled by the momen- 
tum conservation. 

In non-central heavy-ion collisions, the transverse hy- 
drodynamic expansion after the initial thermalization 
leads to a strong azimuthal asymmetry of hadron distri- 



Figure 1: The transverse momentum correlation Ct(<I>i,4>2) 
from the HIJING Monte Carlo model in Au + Au central 
collisions at ^J~s — 200 GeV with (solid) and without (dashed) 
jet quenching. In central collisions, CV(0i,02) only depends 
on ip = <j>\ — cj>2 . The transverse momentum fluctuation in the 
same bin at <p = is subtracted from the correlation. 
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bution or elliptic flow 0, [23| induced by the asymmetry of 
the overlapping geometry and therefore will have similar 
two-hadron correlations as that in Fig. [T] due to mini- 
jets. But such correlation is for the averaged multiplicity 
and with respect to the nuclear reaction plane. This is 
different from the correlation in Eq. ([TJ, which focuses 
on correlation of fluctuation from event to event. Never- 
theless, there might exist contributions from elliptic flow 
to the transverse momentum correlation of our interest. 
Selection of central collisions will therefore help minimize 
the influence of the elliptic flow, which can be further re- 
duced by averaging the transverse momentum correlation 
Ct (01,02) over all possible pairs of bins with fixed width 
80. In central collisions, Ct only depends on ip — (f>\ — (f>2, 
so we can use the notation CV(y) = Ct{4>i ~ <f>2)- 

3. Diffusion equation. After the initial (or partial) 
thermalization, we can assume that the system will obey 
the Bjorken scaling and can be described by the cylin- 
drical coordinates = [r, 77, r, (f>] with metrics g^ v = 
diag(— 1, t 2 , 1, r 2 ), where r = y/t 2 — z 2 is the proper 
time, r\ the spatial rapidity, r the transverse radius and 
(f> the azimuthal angle in transverse plane. We will only 
consider the central rapidity region in central collisions 
and therefore the flow velocity can be cast in the sim- 
ple form = 7t[1, 0, vt,0], where 77- = 1 / \J 1 — v 2 -, and 
vt = vt (t, r, (f>) is independent of the spatial rapidity 77 in 
the Bjorken scaling scenario. Under this cylindrical sys- 
tem, we can neglect the diagonal components of the shear 
stress tensor 7r rr , ir TT and tt^ as compared to their ideal 
counterpart Tg r , Tq t and Tq , and consider the only 
off-diagonal component T^ r = Ti^ r = —r/s-f^d^VT/r 2 of 
the shear stress tensor in the r-component of the fluid 
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0, which can be written as 
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+ i( e + P) 7 |4 + |- [(e + P) 7 |4] , (2) 

where ?y s denotes the shear viscosity and should be inde- 
pendent of the azimuthal angle (f> in the case of central 
heavy- ion collisions, and e and P are energy density and 
pressure, respectively. 

As shown by the HIJING results in Fig. [TJ the mag- 
nitude of the transverse momentum correlation due to 
mini-jets is quite small. Therefore, in the study of the 
evolution of the transverse momentum correlation, the 
fluctuation of energy-momentum density as a result of 
the correlation can be considered as a perturbation in 
a near-equilibrium system. Such a perturbation can be 
characterized by the perturbation 5vt in the transverse 
flow velocity vt —* vt + 5vt- Neglecting fluctuation in e 
and P and considering central collisions where vt is inde- 
pendent of 4>, the equation for the corresponding fluctu- 
ation in the transverse flow velocity can be derived from 
Eq. © as, 
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The above is a diffusion-like equation, which can be 
seen by observing that in the co-moving frame where 
vt = 0, the above equation becomes the standard dif- 
fusion equation 8(tSvt) / dr — D8 2 (t5vt) /d(j> 2 with the 
diffusion constant D — r/ s /(Tsr 2 ) if one neglects the 
time-dependence of the enthalpy density sT. 

With a given initial distribution of 5vt one can solve 
the above diffusion equation and then calculate the final 
transverse energy fluctuation at the freeze-out surface, 
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= Ar) J (v T - jTSgrrrdrdc/), (4) 

where A77 is the spatial rapidity interval, 

dr 

eEu = (-1, 0, —,0)rTdrd(j)dr] (5) 
or 

is the freeze-out hyper-surface area element, and 

Sg T (T, r, 4>) = 7! (e + P)Sv T (6) 

is the fluctuation in transverse momentum density. One 
can then compute the azimuthal correlation of the trans- 
verse energy at freeze-out, 

(^yw r ,r , 

Ct((P) = ,„ 7-7v\2 — / r\T\dr\ / r 2 T 2 dr 2 

V^TW/Out J OUt J OUt 

x {6gT(Ti(rt),r l ,4>)6g T (T2{r2),r2,(l> + <p)) mt , (7) 



where the average should be over an ensemble of initial 
conditions for the spatial distribution of Svt and 
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= AriSifi f ( 1 — vt-^t- \ 7t£ rrdr. (8) 
J out \ or J 

is the averaged transverse energy per bin. 

To solve the diffusion equation ([3]), one has to provide 
the initial condition for Svt as well as the space-time 
evolution of vt, £ and P as functions of r and r. To ob- 
tain the space-time dependence of vt, e and P, we solve 
the ideal hydrodynamic equations in (2+l)-dimension 
numerically [24l [2a | with the initial transverse velocity 
vt{to) = and initial temperature of T = 359 MeV 
(corresponding to energy density eo = 30.0 GeV/fm 3 ) at 
the center of the overlap region of central Au + Au colli- 
sions at an initial time To = 0.6 fm/c. The initial trans- 
verse distributions of e and P are given by the Glauber 
model for nuclear collisions with the Woods-Saxon nu- 
clear distribution. With such initial conditions for (2+1)- 
d ideal hydrodynamic equations, the final transverse en- 
ergy dEr/dri — 667 GeV, calculated at a freeze-out 
temperature T/ = 128 MeV (or energy density 0.075 
GeV/fm 3 ) in central (impact-parameter 6 = 0) Au + Au 
collisions at y/s = 200 GeV/n, is consistent with the 
experimental value rfPy/dr? = 606 ± 35 GeV for 0-5% 
central Au + Au collisions [26| , which is also close to the 
value of HIJING result dE T /dri = 750 GeV. Note that 
the transverse energy as defined in Eq. (1), Eqs. ((1])-© 
and © is slightly different from the calorimetric energy 
measured in experiments. However, the differences are 
small if one considers only hadrons in the central rapid- 
ity region. 

We will use the HIJING model to estimate the initial 
condition for Svt- We use the event ensemble for 5Et 
generated by HIJING (with jet quenching) to extract the 
angular distribution for 5vt(tq.) at the initial time accord- 
ing to Eqs. (0| and ([5]), assuming 5vt(tq) is independent 
of r in the overlapped region for simplicity. This proce- 
dure of extracting the velocity fluctuation from HIJING, 
however, must rely on our knowledge of the initial val- 
ues of energy density e and pressure P at the start of 
the hydrodynamic evolution when the system is in com- 
plete or partial thermalization at To. The initial values of 
energy density and pressure that we use for solving the 
ideal hydrodynamic equation, however, would give twice 
the total transverse energy from HIJING if we assume a 
free-streaming freeze-out at the initial time tq. It is the 
PV work by the Bjorken longitudinal expansion in the 
(2+l)-d ideal hydrodynamics that causes the total trans- 
verse energy at the final freeze-out to be reduced by half. 
To match this initial condition for the ideal hydrodynam- 
ics, we will then multiply 8Et from the HIJING events 
by a factor of 2 so that the transverse momentum corre- 
lation evaluated from Eqs. (HJ)-© with a free-streaming 
freeze-out at the initial time tq becomes identical to the 
HIJING result (with jet quenching) as shown in Fig. 1. 
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Figure 2: (color online) The transverse momentum correlation 
evolved from the diffusion equation as a function of azimuthal 
angle difference tp with different values of n a /s on the freeze- 
out hyper-surface. The EOS is that with the first order phase 
transition or EOSQ. 
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In the calculation of the transverse momentum corre- 
lation Cr(v) from HIJING events as shown in Fig. [TJ 
we have removed the contribution from transverse mo- 
mentum fluctuation Y^A^Ti) ~ (En) 2 in the same bin 
at ip — 0. Such fluctuation should also exist in the ini- 
tial condition for 8vt from HIJING events and will in- 
fluence the transverse energy correlation at (p ^ at a 
later time due to the diffusion according to Eq. We 
will use the following procedure to remove the influence 
of the same-bin transverse momentum fluctuation on the 
final transverse momentum correlation through the dif- 
fusion equation. We first generate an event ensemble of 
Svt whose distribution in is random and the initial 
transverse momentum correlations all vanish except at 
the first bin (tp = 0). The magnitude of the random 
Svt distribution is adjusted such that the initial trans- 
verse correlation in the first bin at tp = matches that 
given by the HIJING events. We then solve the diffusion 
equation ([3]) with such a random distribution of Svt and 
calculate the event-averaged transverse momentum cor- 
relation C^p (ip) at the same final freeze-out surface. Such 
correlation will be subtracted from the results with the 
HIJING initial condition to give us the final transverse 
energy correlation without the influence of the same-bin 
fluctuation. 

Shown in Fig. [2] are the final event-averaged transverse 
momentum correlations Ct(<P>) from the solution of the 
diffusion equation ^ with different values of the shear 
viscosity to entropy density ratio rj s /s. The azimuthal 
dependence of Sgx from the solution of Eq. ([3]) will be 
independent of time for zero viscosity rj s — 0. However, 
its normalization will change because of the longitudinal 
and transverse expansion. During the evolution, both 
the transverse momentum fluctuation and the averaged 
transverse energy will decrease with time. However, the 
latter decreases faster as governed by the ideal hydrody- 
namics. This leads to the increase of the magnitude of the 



transverse correlation Ct(<P>) due to hydrodynamic ex- 
pansion even for ideal fluid with rj s /s = while the shape 
of the correlation remains the same. This is consistent 
with the observation that the magnitude of the trans- 
verse energy correlation from HIJING with jet quenching 
is larger than that without as shown in Fig. 1. 

As one can see, the transverse momentum correlation 
becomes diffused as one increases the value of shear vis- 
cosity to entropy density ratio, which characterizes the 
momentum exchange among nearby fluid elements in dif- 
ferent bins of azimuthal angle. This is consistent with the 
HIJING results with and without jet quenching since jet 
quenching happens at the earliest stage of heavy-ion col- 
lisions during the initial thermalization process where the 
shear viscosity is the largest. The diffusion of the correla- 
tion in the near-side due to the viscous evolution is much 
stronger than the away-side. We have also studied the 
dependence of the final transverse momentum correlation 
on the equation of state (EOS). We find that the shape 
of the final correlation is not very sensitive to the EOS, 
so we present in Fig. [5] only the result of the EOS with 
the first order phase transition or EOSQ. 

4- Conclusions and discussions. We have studied 
transverse energy correlation in azimuthal angle and its 
evolution in high-energy heavy-ion collisions. We show 
that such correlation will have a characteristic structure 
of two Gaussian peaks in high-energy heavy-ion collisions 
at the RHIC energy and above due to the dominance of 
mini-jets in the initial parton production. We find that 
within the HIJING Monte Carlo model, the transverse 
energy correlation will be significantly broadened or dif- 
fused by jet quenching in the early stage of thermalization 
in heavy-ion collisions. Assuming that the residual cor- 
relation still exists in the nearly-equilibrated system, we 
derived a diffusion-like equation for the local transverse 
energy density fluctuation by taking linearized approxi- 
mation of the viscous fluid equations. With the initial 
condition as given by the HIJING results, we solved such 
a diffusion equation in combination with the space-time 
evolution of the averaged thermodynamic variables and 
flow velocity as given by the ideal hydrodynamic equa- 
tions. The final transverse energy correlation of produced 
hadrons averaged over many events is found to be diffused 
because of the shear viscosity. For larger values of the 
shear viscosity to entropy density ratio, r] s /s ~ 0.2 — 0.4, 
the diffusion is as significant as that caused by jet quench- 
ing. One can therefore essentially study the degree of 
mini-jet thermalization in heavy-ion collisions as well as 
the viscous effect of the strongly coupled QGP [U [H[ 
by measuring the transverse momentum correlation in 
heavy- ion collisions in comparison with that in p + p and 
P + A. 

The correlations of transverse momentum per parti- 
cle in both pseudo-rapidity and azimuthal angle have 
been measured by the STAR experiment [2!| at RHIC 
through an autocorrelation technique which inverts the 
bin-size dependence of the event-wise mean transverse 
momentum fluctuation to extract the bin-bin correlation. 
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Such correlation is different from the correlation of the 
transverse momentum which is the total transverse mo- 
mentum within a bin with given size S(f>. Nevertheless 
the transverse momentum correlation that is inferred in 
the STAR experiments and its centrality dependence [29[ 
might point to medium modification of the jet fragmenta- 
tion and hadronization process that could also influence 
the transverse momentum correlation as studied in this 
paper. Such effect beyond the normal jet hadronization 
and hydrodynamic evolution of the medium will be inter- 
esting to investigate through experimental measurement 
of the transverse momentum correlation as proposed in 
this paper. 
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